Quantitative modeling, optimization, and verification of 
°3Ni-powered betavoltaic cells based on three-dimensional 


ZnO nanorod arrays 

Zan Ding!:”, Tong-Xin Jiang!:*, Ren-Rong Zheng! ?, Na Wang! ?, Li-Feng Zhang’, Shi-Chao 
Liu’, Xin Li*: **, Hai-Sheng San (https://orcid.org/0000-0002-6019-923X)!:* * 

'Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, 
Xiamen 361005, China 

2 Shenzhen Research Institute of Xiamen University, Shenzhen 518000, China 

3 China Institute of Atomic Energy, Beijing 102413, China 

* Shanghai Institute of Space Power-Sources, Shanghai 200245, China 

= ** Corresponding authors, sanhs @xmu.edu.cn (H. S. San); lixin_0128 @sina.com (X. Li) 
This study was supported by the National Natural Science Foundation of China (No. 
12175190), the Special Funds for Central Government Guiding Shenzhen Development in 


Science and Technology (No. 2021Szvup066) and the Natural Science Foundation of Fujian 
Province, China (No.2022J02006). 


Abstract 


Betavoltaic cells (BCs) are promising self-generating power cells with long life and high 
power density. However, the low energy conversion efficiency (ECE) has limitations in 
practical engineering applications. Wide-bandgap semiconductors (WBGSs) with 
three-dimensional (3-D) nanostructures are ideal candidates for increasing the ECE of BCs. 
This paper proposes hydrothermally-grown ZnO nanorod arrays (ZNRAs) for ®Ni-powered 
BCs. A quantitative model was established for simulation using the parameter values of the 
dark characteristics, which were obtained from the experimental measurements for a 
simulated BC based on a Ni-incorporated ZNRAs structure. Monte Carlo (MC) modeling and 
simulation were conducted to obtain the values of the f energy deposited in ZNRAs with 
different nanorod spacings and heights. Through the simulation and optimization of the 3-D 
ZNRAs and 2-D ZnO bulk structures, the performance of the “Ni-powered BCs based on 


both structures was evaluated using a quantitative model. The BCs based on the 3-D ZNRAs 


structure and 2-D ZnO bulk structure achieved a maximum ECE of 10.1% and 4.69%, 
respectively, which indicates the significant superiority of 3-D nanostructured WBGSs in 


terms of increasing the ECE of BCs. 
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1. Introduction 


With the rapid development of wireless sensing networks for harsh environments (deep 
sea, deep ground, deep space, and polar and desert regions), the demand for micropower 
sources has increased significantly. Compared with conventional cells (chemical cells, fuel 
cells, and solar cells), betavoltaic cells (BCs) with longer lifetimes, stronger environment 
adaptiveness, and higher energy density have good application prospects and remarkable 
market potential [1]. 

Recently, reviews on BCs were published by Olsen et al. [2], Prelas et al. [3], and 
Spencer et al. [4]. According to Olsen’s calculation model, the possible efficiency cap can 
reach as high as 35% for BCs employing wide-bandgap semiconductors (WBGSs). However, 
the recorded energy conversion efficiency (ECE) reported for experimental betavoltaics is far 
less than 35% when using the available WBGS. In previous decades, research on BCs has 
received more attention, and related experimental studies have recorded the ECE in the order 
of 6% when using available WBGSs (4H-SiC and diamond) [5-9]. To increase the output 
power (Pout) of BCs, multiple BCs are stacked modularly to generate microwatt-level output 
power [10]. However, the issue of low ECE has not been resolved to facilitate the practical 
engineering applications of BCs. It is expected that the loading and utilization efficiency of 
radioactive sources in energy conversion structures (ECSs) can be further improved using 
new strategies. 

It is well understood that an ECS with a larger specific surface area can load more 
radioactive sources, which means that BCs can generate higher output power. In conventional 
planar diode structures, the ECE and output power density are limited owing to the low active 
area in the planar structure. In contrast, three-dimensional (3-D) nanostructures with a larger 
specific surface area can significantly increase the J energy conversion, which can be 
explained by the following four basic facts: i) the filling of radioisotope material into the 
interspace of 3-D nanostructures results in significant increase in the loading amount of 
radioisotope materials; 11) the combination of a radioisotope source with 3-D nanostructures 
results in the interaction of p particles with nanostructured ECSs in all directions, and thus the 


collection efficiency of the p particles dramatically increases; iii) the self-shielding effects of 


the radioisotope source are cancelled owing to the nanometer-level thickness of the 
radioactive structure; iv) the effective junction region area of the 3-D nanostructure is 
significantly extended. 

Recently, several evaluation studies considering a 3-D nanostructure have been 
conducted based on theoretical analysis or experimental verification. For example, Sun et al. 
[11] developed a BC based on a 3-D porous silicon p-n diode structure loaded with tritium gas, 
and achieved an ECE of 0.22%, which was ten times larger compared with that of a 2-D 
planar diode structure. Murphy et al. [12] designed a BC based on 3-D silicon diode structure 
incorporated with '47Pm source, which achieved an ECE of 5.8% and power density of 25 
mW/cm*. In our previous studies, WBGS TiO: (bandgap = 3.2 eV) and ZnO (bandgap = 3.37 
eV) nanostructures were used to fabricate the ECSs of BCs [13-19]. A cost-effective anodic 
oxidation method and a hydrothermal method were used to prepare highly ordered and 
free-standing nanotube and nanorod array structures, respectively. Using a Nickel-63 (Ni) 
sheet as the radioisotope source and electrode, an enhanced betavoltaic effect was observed in 
the BCs. In the last decade, 3-D ZnO nanorod arrays (ZNRAs) have attracted extensive 
attention owing to their large bandgaps, high aspect ratios, and large specific surface areas. 
Compared with TiO2, ZnO has longer electron lifetime and higher electron mobility [20], 
which is considered as one of the most promising electrode materials for photoelectronic 
devices. For example, in photodetectors using ZNRAs as the photoelectrodes, electrons can 
be rapidly transported along ZnO nanorods, which helps in suppressing the recombination of 
electron-hole pairs (EHPs) [21]. Additionally, ZNRAs can be synthesized on a wide variety of 
substrates for versatile electrode preparation [22]. However, devices based on pure ZNRAs 
are subjected to high carrier recombination owing to the lack of an effective separation 
mechanism for the EHPs. For BCs using ZNRAs, an effective method is to incorporate a 
metal-based radioactive material, such as “Ni or “Pm, into the interspaces of ZNRAs to 
form a core-shell structure with a Schottky junction that can considerably separate the EHPs. 
Additionally, the loading amount of the radiation source can be increased significantly 
compared with using Ni sheets as the radiation source, which is conducive to increasing the 
Pout of BCs. 


Recent studies on the optimal design of BCs with 2-D bulk structures have shown that 
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the size and structure of ECSs and the thickness of the radioactive source have significant 
influence on the performance of BCs [23-29]. However, few studies have been conducted on 
the influence of 3-D nanostructured ECSs combined with radioactive sources on the 
performance of BCs. Therefore, an accurate theoretical model must be established to evaluate 
the performance of BCs based on 3-D nanostructured WBGSs. 

This paper proposes a theoretical model for evaluating the performance of BCs based on 
3-D nanostructured WBGSs. A freestanding ZNRAs structure combined with normal Ni was 
used to simulate the ECSs to obtain the parameters of the dark characteristics. A Monte Carlo 
(MC) model based on the MCNP code was established to investigate the interaction between 
incident 2 particles and the ECSs. By combining the theoretical model with experimental data, 
the BCs based on the “Ni-incorporated ZNRAs structure can be accurately evaluated and 


further optimized. 


2. Theory, design and experiment 


2.1. MC simulation for p energy deposition on 3-D nanostructure 


(b) 
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Fig. 1. (Color online) Schematic diagram of BCs based on (a) 2-D ZnO bulk structure and (b) 
3-D ZNRAs structure. 


To design an optimal structure of BCs, MC simulation [30] was conducted to investigate 
the interaction of £ radiation with 3-D “Ni-incorporated ZNRA structures. As the £-emitting 
source, this study selected “Ni with a half-life of approximately 100 years, average energy of 
17.4 keV, and maximum energy of 66 keV. Two types of BCs were designed to compare the 
energy deposition simulations. As shown in Fig. la, a sandwich-type BC based on a 2-D ZnO 
bulk structure consisted of an Al-doped ZnO (AZO)/glass substrate, ZnO film, and “Ni/Ni 


layer deposited on the ZnO film. In contrast, the 3-D nanostructured BC consisted of an 


AZO/glass substrate, 3-D ZNRAs combined with “Ni, and a top Ni electrode (Fig. 1b). 
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Fig. 2. (Color online) MC modeling and simulation of 2-D ZnO bulk structure and 3-D 
ZNRAs structure under £ irradiation. Scattering trajectories of 2 particles in (a) 2-D ZnO bulk 
structure and (e) 3-D ZNRA structure; energy deposition distributions of p particles in (b—c) 
2-D ZnO bulk structure and (f-g) 3-D ZNRA structure on y-z plane and x-y plane; 
dependence of energy deposition density and deposition ratio on depth in (d) 2-D ZnO bulk 


structure and (h) 3-D ZNRAs structure. 


To qualitatively understand the interaction of incident £ particles with ECSs, 2-D ZnO 


bulk and 3-D ZNRAs structures were modeled and simulated using the MC method. In this 


mode, the ZnO thickness was set to 3 um, and the diameter and length of the ZnO nanorods 
were set to 300 nm and 5 um, respectively. An electron beam with a diameter of 1.8 um and 
incident angle of 45° along the horizontal direction was used as the exciting source, and 
consisted of 1000 £ particles with an average energy of 17.4 keV for the particles of “Ni. 
The scattering trajectories and energy deposition distribution of the # particles in both 
structures were simulated using the CASINO software, as shown in Fig. 2. As shown in Fig. 
2a, a fraction of the incident £ particles was reflected from the surface of the bulk ZnO, which 
led to certain energy loss in the 2-D ZnO bulk structure. The energy deposition distributions 
of the incident 2 particles in the 2-D ZnO bulk exhibited an incident depth of approximately 
1.47 um on the y—z plane (Fig. 2b) and an interaction radius of 2.9 um on the x-y plane (Fig. 
2c). Fig. 2d shows the dependence of the energy deposition density and deposition ratio on the 
depth of the 2-D ZnO bulk structure. The maximum depth of energy deposition is 
approximately 1.47 um, and the energy deposition ratio reached 32%. In contrast, the MC 
simulation yielded different results for the ZNRAs. The scattering behaviors of the 2 particles 
in the ZNRAs imply a higher structure trapping for the 2 particles compared with the 2-D 
ZnO bulk structure (Fig. 2e). Moreover, the simulation results for the energy deposition 
distribution in the ZNRAs (Fig. 2f and 2g) reveal that incident energy can be deposited in a 
very deep range along the length of the nanorods owing to the multiple scattering and 
absorption of p particles. This can be further verified by the simulation results for the 
dependence of the energy deposition density and energy deposition ratio on the depth in 3-D 
ZNRAs structures, as shown in Fig. 2h. Compared with the 2-D ZnO structure, the energy 
deposition depth in the 3-D ZNRAs structure could reach approximately 4.5 um with an 
energy deposition rate of approximately 85%. This result indicates that the 3-D ZNRAs can 
harvest more J energy in the ECSs compared with the 2-D ZnO structures, which results in 


the higher ECE and Pou of BCs. 
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Fig. 3. (Color online) (a) Schematic energy band diagram of BC based on °Ni@ZNRAs and 
(b) transferring process of 2 excited carriers in “Ni@ZNRAs; (c) energy deposition and 
deposition ratio of incident £ particles of devices based on “Ni@ZNRAs with different 
nanorod spacings; (d) energy deposition distribution and deposition ratio of incident ~ 


particles of devices based on °Ni@ZNRAs with different nanorod heights. 


Figure 3 shows the working principle of the BCs based on “Ni@ZNRAs. A schematic 
energy band diagram is shown in Fig. 3a. An extended 3-D heterojunction of “Ni/ZnO is 
formed with a Schottky contact (®spu = Wm — y) in the interface of °Ni/ZnO and an ohmic 
contact in the interface of AZO/ZnO. The £ particles emitted from the Ni source are incident 
on the ZnO lattice, where impact ionization results in the energy deposition and generation of 
electron-hole pairs (EHPs). The EHPs are separated by a built-in electric field in the depletion 
region of the “Ni@ZnO Schottky junction. In this process, electrons migrate from the 
conduction band (with energy level Ec) of ZnO to the AZO substrate, while holes are 
transferred from the valence band (with energy level Ev) of ZnO to the “Ni electrode [31]. 


Figure 3b shows the horizontal EHP separation and longitudinal carrier transport. 


To simulate the 8 energy deposition in the ZNRAs with satisfactory accuracy, the MCNP 
software was used to establish the MC model of the °Ni@ZNRAs. This study used “Ni as an 
isotropic volumetric source integrated into the gap of the nanorods, and the number of 
incident 2 particles was set to 1x10°. The 8 energy deposition in the ZNRAs by radiation 
from the “Ni source was simulated using the actual spectrum of the emitted electrons. The 
energy deposition in the ZNRAs was calculated using the F8 pulse counting in the MCNP 
program, and the deposited energy could then be employed to calculate the actual current 
density using Eq. (2), as described in Sect. 2.2. The dependence of the energy deposition and 
deposition ratio of the incident £ particles on the spacing and height of the nanorods is shown 
in Fig. 3c and 3d, respectively. As the spacing of the nanorods increased, the loading amount 
of Ni increased, which led to an increase in the total -emitting energy. However, the 
number of nanorods per unit area decreased as the spacing increased. Consequently, the 
energy deposition rate decreased. The trade-off reached a balance at the spacing of 
approximately 200 nm, where the maximum energy deposition was achieved (Fig. 3c). As the 
nanorod height increased, as shown in Fig. 3d, both the deposited energy and energy 
deposition ratio increased, and the energy deposition rate exhibited a trend of slow increase 
when the height exceeded 4 um. 

2.2. Quantitative modeling for BCs based on 3-D ZNRAs 
A quantitative model should be established to accurately evaluate the performance of the 
BCs. Because the working principle of BCs is similar to that of solar cells, a conventional 
theoretical model of solar cells can be improved to establish a model of 3-D nanostructured 
BCs. As high-energy incident 2 particles on the ECSs, multiple electron-hole pairs (EHPs) 
are generated by the impact ionization of the 2 particles with the semiconductors. The 
EHPs are effectively separated and drawn as the theoretical maximum current density (/max), 


which can be calculated as follows [32]: 


Tax =N g(Ey/V) (1) 
y =2.8E,+ 0.5 eV (2) 


where wy is the ionization energy that generates an EHP; Eg is the bandgap of the 


semiconductor; Nz is the -emitting flux (Ng = A x 3.7 x 10’); A is the activity (mCi) of the 
radioactive source; Ez is the average energy of the 8 particles. The radiation-generated current 


(Ig) in the BCs is given as follows: 


L, = 90-1 p)1-) Oy, = OLE, 1y 3) 


where y is the reflection coefficient of the 2 particles from the semiconductor surface; u is the 
interaction coefficient between the external medium (air or electrode) and incident £ particles; 
a is the self-absorption coefficient of the radioactive source, which is related to the thickness 
of the radioactive source; gis the scattering angle coefficient of the radioactive source, which 
is the ratio of the number of incident particles on the semiconductor to the total number of 
emitted £ particles. For 3-D ZNRAs, y= 0, u= 0, m= 1, e is the electron charge (e = 1.6 10°!” 
C), and EF; is the 2 energy deposited in the semiconductors; Q is the collection efficiency, 
which is the ratio between the number of effective carriers for forming a current in the devices 
and the total number of carriers produced by the impact ionization of the incident 2 particles 
with the semiconductors. Here, Q is the collection efficiency of carriers in the space-charge 
region of the Schottky junction, which is approximately 100%; the collection efficiency of the 


carriers outside of the space charge region can be calculated as follows [33]: 
Q = 1-tanh(d/L) : (4) 


where d is the distance of the carrier moving to the boundary of the depletion region and L is 
the minority carrier diffusion length. 


The width of the Schottky junction space-charge region can be calculated as follows [34]: 


W = ./28,€.(®,,, -V)/(eN,) (5) 


where sois the dielectric constant of the vacuum; gs is the relative dielectric constant of the 
semiconductor; Np is the ZnO carrier concentration (approximately 2.1x10!5 cm™); ®spu is 
the height of the Schottky barrier of Ni@ZnO (ssu = Wm — vy = 0.8 eV), which is, the 
theoretical difference between the work function of Ni (Wm = 5.15 eV) and the electron 
affinity energy of ZnO (y = 4.35 eV). 

From Eq. (5), the width of the depletion region of the Ni@ZnO Schottky barrier can be 


calculated as approximately 6.0410"! um. Because the average diameter of the ZnO nanorods 


(approximately 0.3 um) grown using the hydrothermal method is smaller than the depletion 

region width, Q is approximately equal to 1 for BCs based on the °Ni@ZNRAs structure. 
Figure S1 shows the equivalent circuit model of the BCs (Supporting Information). The 

current density (J) delivered to the load resistor (RL) at voltage V can be calculated as follows 


[35]: 


I =I, —6I,lexp(eV + IaR,)/ (nKT)) -1]-(V +1aR,)/ Ry, (6) 


5=S/S,, a=(D,H)\\(DH,), (7) 


where Rs and Rsn are the series and shunt resistances of the equivalent circuit of the BC, 
respectively. For the specific BC based on 3-D ZNRAs that was prepared for the subsequent 
experiment, Jo, n, and Rs were extracted from its J-V curve measured under dark conditions, 
and the ideal Rsn was set to infinity. The greater height and smaller density of the nanorods 
will result in greater Rs. Considering that both J, and Rs are variables related to the structural 
parameters of 3-D ZNRAs, such as the area of the Schottky junction (S), nanorod array 
density (D), and nanorod height (H), two factors (ô and a) are combined with /, and Rs in Eqs. 
(6-7). Thus, Ho and Ds are the structural parameters of the specific 3-D ZNRAs structure. 

The maximum short-circuit current (Zs) of the BCs was calculated using Eq. (6). The 


open circuit voltage (Voc) of the BCs can be calculated as follows: 


V =nkT in, /1,+1)/e, (8) 


where k is the Boltzmann constant (k =1.38 x 10° J/K), T is the thermodynamic temperature 
at room temperature (T = 300 K), Jo is the reverse saturation current, and n is an ideal factor. 
The maximum output power (Pmax), filling factor (FF), and energy conversion efficiency 


(n) can be expressed as follows [36]: 


P „=I VFF 
max sc" oc (9) 

FF = (1 Vn) UVa) (10) 

n = Px | P, x100% = P x / (eN, E,) x 100% (11) 


where Jmp and Vmp are the current and voltage, respectively, corresponding to Pmax; Pg is the 


total power of the radioactive source; No is the total amount of £8 particles (No = 3.7x10’x@, 


where @ is the activity of the radioactive source) emitted from the radioactive source. 


2.3. Preparation of BCs based on 3-D ZNRAs nanostructure 

The preparation process of the Ni-incorporated ZNRAs structure is shown in Fig. 4a. The 
ZNRAs were grown on Al-doped ZnO (AZO) conductive glass (Ro< 6 Q, Nanotech, China) 
using a hydrothermal method [16]. First, the AZO conductive glass was cleaned using 
CH3COCH3, ethanol (C2HsOH), and deionized water in an ultrasonic bath for 15 min. After 
drying, the specified electrode region in the AZO glass substrates was covered with high 
temperature tape, and the substrates were then transferred to a reactor with 10 ml of growth 
solution consisting of 0.05 M Zinc nitrate hexahydrate [Zn(NO3)2°6H20], 0.05 M 
hexamethylenetetramine (CH2)6N4, HMTA), 0.02 M polyetherimide (PEI), and deionized 
water. The temperature of the hydrothermal treatment was maintained at 95 °C for 4 h. Next, 
the samples were cleaned with deionized water and dried at 150 °C for 10 min. Finally, the 
samples were soaked in 10% concentration of hydrogen peroxide (H202) for 1 min for surface 
passivation, and then annealed in air at 400 °C for 2 h to obtain the ZNRAs [37]. Passivation 
treatment with H2O2 can reduce the surface defects of the ZnO nanorods, which favors the 
improvement of the Schottky diode characteristics. The incorporation of Ni into the ZNRAs 
was performed using a direct current (DC) magnetron sputtering technique. The vacuum 
chamber was evacuated to a pressure of 5x10% Pa, and Ar gas was then introduced into the 
chamber to reach the constant pressure of 0.6 Pa. Subsequently, Ni nanoparticles were 
deposited onto the ZNRAs using a DC power of 400 W at room temperature for 9 min. Fig. 
4b shows a 3-D explosive diagram of the device based on the Ni-incorporated ZNRAs 
structure. To perform the /-V measurements, a planar Ni sheet was mechanically compressed 
onto the surface of the as-prepared samples to form a sandwich-type device with a 


Ni/Ni@ZNRAs/AZO/glass structure. 
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Fig. 4 (Color online) (a) Schematic illustration of preparation process flow of 
Ni-incorporated ZNRAs structure; (b) 3-D explosive diagram of device based on 
Ni-incorporated ZNRAs structure; (c) photographs of (1) AZO glass, (11) ZNRAs/AZO glass, 


and (iii) Ni-incorporated ZNRAs/AZO/glass samples. 


2.4. Material characterization and device measurements 

Scanning electron microscopy (SEM, ZEISS microscope, Germany) and 
energy-dispersive X-ray spectroscopy (EDS) were used to analyze the morphology and 
chemical composition of the samples. The crystal structures were characterized by X-ray 
diffraction (XRD, D8 Advance, Bruker, Germany), and 7- V measurements were performed in 


a Faraday cage using a Keithley Model 2450 source meter. 


3. Results and discussion 


3.1. Experimental results and discussion 
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Fig. 5. (Color online) (a) XRD spectra of Ni@ZNRAs structure; (b) cross-section-view and (c) 
top-view SEM images of ZNRAs; (d) cross-section-view and (e) top-view SEM images of 
Ni@ZNRAs structure; (f-h) EDS mappings of Zn, Ni, and O. 

The XRD pattern of the Ni@ZNRA structure is shown in Fig. 5a. All diffraction peaks 
can be indexed to the hexagonal structure of wurtzite ZnO (JCPDS card No. 36-1451) and 
cubic structure of Ni (JCPDS card No. 70-1849). The typical 20 peaks at 34.4° and 36.3° are 
assigned to the (002) and (101) planes of ZnO, respectively. A new diffraction peak was not 
observed in the Ni@ZNRA structure, which indicates that the introduction of Ni did not 
change the crystal structure of the ZnO nanorods. The growth direction of the (002) crystal 
plane of ZnO is preferred to the c-axis orientation perpendicular to the substrate, as can be 
verified from the cross-sectional and top-view SEM images of the ZNRAs (Fig. 5b and Fig. 
5c). The average diameter and height of the nanorods is 300 nm and 4 um, respectively. The 
cross-section-view and top-view SEM images of the Ni@ZNRAs are shown in Fig. 5d and 5e, 
respectively. As can be clearly observed, Ni particles were uniformly filled into the gap 
between the ZnO nanorods and stacked to form a continuous Ni layer covering the surface of 
the ZNRAs. The SEM-EDS mappings show that Ni was adequately distributed on the surface 


of the ZnO nanorods, forming a core-shell structure (Fig. 5f-h). 


( 


) 1.5x10* 


(b) 10* 


10° 


(c) 


wE -t fo=4.2x10' Alem? 


dV/d(Inl) 


Current density (A/cm?) & 
P D y 
A 
C] 
Current density (A/cm?) 


-3.0x10° B 0.0 
-25 -2.0 -1.5 -1.0 05 0.0 05 10 15 20 25 0.0 0.5 1.0 1.5 2.0 25 3.0 0.0 5,0x10° 1.0x10° 1.5x10° 2.0x10° 2.5x10* 3.0x10° 


Voltage (V) Voltage (V) Current density( A/cm?) 


Fig. 6. Dark characteristics of devices based on Ni@ZNRAs/AZO structure: (a) linear-scale 
I-V curve in bias voltage range of -2.5—2.5 V; (b) logarithm-scale /-V curve in forward bias; (c) 
dV/d(In/)-/ curve. 


Figure 6 shows the dark characteristics of the devices based on the Ni@ZNRAs/AZO 
structure. A typical J-V curve of the Schottky diode is shown in Fig. 6a, confirming the 
existence of the Ni/ZnO Schottky junction. The reverse saturation current of the device is 


4.2x10!! A/cm”, as determined by extrapolating the linear region of the logarithm-scale /-V 


curve in the forward bias to zero bias voltage (Fig. 6b). The ideal factor n and series resistance 
Rs of the Schottky diode can be extracted from the dependence of dV/d(InJ) on the current (Fig. 
6c) [38]. The extracted n and Rs are 6.96 and 248.45 Q, respectively. This ideal factor n is 
larger than 2, which indicates that the junction quality of Ni@ZNRAs is not perfect owing to 
the crystal defects or insufficient filling of Ni in the ZNRAs [39-40]. Fig. S2 (Supporting 
Information) shows the J-V curve of the devices based on the Ni@ZNRAs/AZO structure 
before H202 passivation and the calculated n value (n = 11.2). As can be seen, the n value 
before passivation is significantly larger compared with that after passivation, indicating that 
H20>2 passivation can significantly improve the diode quality. The presence of a large amount 
of interface defects in the ZnO nanorods may deteriorate the rectification characteristics of 
Schottky devices, leading to an increase in the reverse saturation current and ideal factor. In 
turn, the increase in reverse saturation current will decrease the open-circuit voltage, resulting 
in the decrease of output power, which can be calculated using Eq. 8. This indicates that H202 
passivation improves the output power of the BCs. To evaluate the electrical performance of 
BCs based on the “Ni@ZNRAs structure, the above-mentioned parameter values of the 


Schottky diode are used in the quantitative model of BCs, as described in Sect. 2.2. 


Table 1. Parameter values used in quantitative model of BCs. 


Parameters Fixed value Variable value 
Reverse saturation current, Io (A/cm?) 4.2x1071!! — 
Ideal factor, n 6.96 — 
Series resistance, Rs (Q) 248.45 — 
Nanorod diameter, D (nm) 300 — 
Nanorod height, H (um) — 0.5—10 

Width of nanorod gap, W (nm) — 40—400 
Specific activity of Ni, Asp (Ci/g) 57 — 
Average energy of “Ni, Es (keV) 17.4 — 
Mass density of “Ni, p (g/cm?) 8.902 — 


Table 1 lists the parameter values applied to the quantitative model of BCs. Notably, 


these parameter values can be controlled and adjusted when growing ZNRAs using the 
hydrothermal method. Tao et al. [41] developed a combined lithography-assisted and 
hydrothermal growth approach to selectively control the spacing and size of ZNRAs on the 
ZnO seed layer. Sinem et al. [42] prepared ZNRAs with nanorod lengths greater than 60 um 
by adjusting the contents of ammonium hydroxide and polyethyleneimine using a 
hydrothermal synthesis method. This means that the optimal values of nanorod parameters, 
such as the diameter, height, and density, can be prepared by adjusting process parameters 
such as the density of the ZnO seeds, concentration of the reactants, time and temperature of 


reactions, and pH values [43]. 


3.2. Simulated results and discussion 

Using the fixed parameter values listed in Table 1 and an average nanorod height of 4 
um, the calculated J-V and P-V characteristics of BCs based on the °Ni@ ZNRAs structure 
with different nanorod spacings are as shown in Fig. 7a and Fig. 7b, respectively. By 
extracting Pmax from Fig. 7b and calculating the ECE according to the total activity of “Ni 
incorporated into the gap of the ZNRAs, the dependences of Pmax and ECE on the variable 
spacing of the nanorods can be obtained as shown in Fig. 7c. As can be seen, as the nanorod 
spacing increases, the activity of the “Ni incorporated into the gap of the nanorods also 
increases, resulting in the rapid increase of output power. However, as the nanorod spacing 
increases further, the nanorod density is reduced, which reduces the energy deposited on the 
ZNRAs and thus reduces the output power. Moreover, a maximum Pou at 875.4 nW/cm? 
corresponding to Jsc=676.2 nA/cm? and Voc=1.85 V was achieved with FF=0.699 and 
n=6.07%. In contrast, Figure 7c shows a monotonic decline in the ECE curve as the nanorod 


spacing increased. The ECE reached a maximum value of 9.3% when the spacing was set to 
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Fig. 7 (Color online) (a) /-V and (b) P-V characteristics of devices based on @Ni@ZNRAs 
with different nanorod spacings; (c) dependences of Pmax and ECE of devices based on 
°3Ni@ZNRAs on nanorod spacing; (d) Z-V and (e) P-V characteristics of devices based on 
SNI@ZNRAs with different nanorod height; (f) dependences of Pmax and ECE of devices 
based on °Ni@ZNRAs on nanorod height. 


The effects of the nanorod height on the ECE and Pmax of the BCs were investigated under 
the optimal nanorod spacing. Fig. 7d and 7e shows the calculated /-V and P-V characteristics 
of the BCs based on the “Ni@ ZNRAs structure with different nanorod heights. As can be 
seen, /sc increased with the nanorod height, while Voc was kept at 1.62 V without any change. 
Fig. 7f shows the dependence of the Pmax and ECE of the BCs on the nanorod height. The 
energy deposition ratio increased slowly when the nanorod height exceeded 4 um (Fig. 3d), 
resulting in a slow increase in the saturation level of the energy conversion efficiency. As the 
nanorod height increased, the ECE tended toward the saturation value of approximately 10%, 
and the maximum Pox of 1185.9 nW/cm? was achieved when the height was 10 um, 
corresponding to a loading amount of 114.37 mCi, Js: of 1081.2 nA/cm’, and FF of 0.673. 
Although long nanorods are beneficial for increasing the output power, it is difficult to fully 
deposit Ni into a high-aspect-ratio ZNRAs structure in practical operation. Hence, 


experimental investigation is required to determine the practical deposition depth of “Ni on 
the ZNRAs. 

To verify the advantages of 3-D nanostructures in enhancing the betavoltaic effect, a BC 
based on a 2-D ZnO bulk structure was investigated for comparison to the BCs based on a 
3-D ZNRAs structure, as shown in Fig. la. According to the MC simulation, the energy 
deposition depth of incident 2 particles in 2-D ZnO bulk is approximately 1.5 um. Moreover, 
the dependence of the energy deposition and the deposition ratio of the incident 2 particles on 


the thickness of “Ni are shown in Fig. 8a. As the thickness of “Ni increases, the 


self-absorption effect of the radiation source shields the 2 particles emitted from the deep 
location of the radiation source, resulting in deposition ratio reduction. As can be seen, the 
energy deposition tended to saturate when the thickness exceeded 1.5 um. The dependences 
of Pmax and ECE of the BCs based on 2-D ZnO bulk on the thickness of “Ni were 
investigated, as shown in Fig. 8b. As the thickness of Ni increased from 50 nm to 150 nm, 
the energy conversion efficiency increased from 4.05% to 4.69% owing to the deposited 
energy increase. As the thickness of “Ni increased further and exceeded 150 nm, the 
self-shielding effect of the radiation source became significant, which reduced the energy 
deposition ratio and thus reduced the energy conversion efficiency. It can be observed that 
Pmax tended toward saturation (approximately 135.2 nW/cm?) when the “Ni thickness reached 
2 um. Beyond this thickness, the self-shielding effect of the radioactive source limits the 
continued increase of output power. Moreover, the maximum ECE reached 4.69% when the 


63Ni thickness was 150 nm, as expressed by Eq. S8 (Supporting Information). 
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Fig. 8. (a) Energy deposition and deposition ratio of incident 2 particles of devices based on 
2-D ZnO bulk with different Ni thickness; (b) dependences of Pmax and ECE of BCs based 
on 2-D ZnO bulk on “Ni thickness. 

The performance comparison between the BCs based on 3-D ZNRAs and those based on 
2-D ZnO bulk is presented in Table 2. The BCs based on the 2-D ZnO film generated higher 
Voc (1.87 V), which is attributed to the generation of lower reverse saturation current 
compared with the 3-D ZNRAs. Owing to the low utilization efficiency of the radioactive 
source, the 8 energy deposited in the semiconductor decreased dramatically, resulting in lower 


Ise (28.1 nA/cm?) and Pmax (36.8 nW/cm’), which led to a low ECE of 4.69% compared with 


the BCs based on the 3-D ZNRAs. According to the Shockley—Quaiser limit of photovoltaics 
[44], the efficiency cap of BCs is predicted to be as high as 13%, as can be calculated by y = 
VəæFF/wy. This means that 3-D ZNRAs with optimal structural design are significantly 
superior and can be used to prepare BCs close to their upper limits. 


Table 2. Comparison of electrical parameters on 3-D and 2-D structured BCs. 


Structures 
Electrical parameters 
3-D ZNRAs 2-D ZnO bulk 
Maximum short circuit current density 
1081.2 28.1 
(nA/cm?) 
Open circuit voltage (V) 1.62 1.87 
Filling factor 0.673 0.701 
Maximum output power (nW/cm?) 1185.9 36.8 
Energy conversion efficiency (%) 10.1 4.69 


4. Conclusion 


In summary, 3-D nanostructured WBGSs are ideal candidates for increasing the ECE and 
output power of BCs. To evaluate and verify their advantages as a highly efficient 
f-energy-conversion structure, this study established a quantitative model for °Ni-powered 
BCs based on freestanding ZNRAs, and carried out MC modeling and simulation to obtain 
the deposited energy and maximum energy deposition depth. In the experiment, a BC based 
on a Ni-incorporated ZNRAs structure grown using the hydrothermal method was simulated 
to obtain the parameter values of the dark characteristics. By optimizing and simulating the 
BCs based on the 3-D ZNRAs structure and 2-D ZnO bulk structure, the performance of the 
°3Ni-powered BCs based on both structures was evaluated using the quantitative model and 
experimental parameter values, respectively. In this study, the optimal nanorod spacing of the 
BCs based on 3-D ZNRAs was 40 nm and 200 nm for maximum energy conversion efficiency 
and maximum output power, respectively. The BCs based on 3-D ZNRAs achieved a 


maximum ECE of 10.1% and maximum output power of 1185.9 nW/cm’, respectively. In 


contrast, the values for the BCs based on the 2-D ZnO bulk are only 4.69% and approximately 
135.2 nW/cm’, respectively. The results confirm the significant superiority of the 3-D 


nanostructured WBGSs in terms of increasing the ECE and output power of BCs. 
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